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In the aftermath of the disappearance of the Malaysian 370 (MH370) flight in March 2014, 
new positioning methods were employed to establish the search area. In the absence of all 
other positioning technologies (Transponder, Radio communications, Radar), these innova- 
tive methods are based on the handshake signals between an INMARSAT satellite and the 
satellite transceiver on board the aircraft. The log of these signals was made public in order 
for the scientific community to engage in solving the mystery of the MH370 trajectory. The 
log indicates the delay between the interrogation and response signals, as well as the relative 
velocity indications, based on the shift of the carrier frequency due to the Doppler-Fizeau 
effect. This paper puts forward an original, independent and accurate positioning method 
and allows the calculation of the MH370 trajectory considering the wind vector field that 
day, the accurate satellite orbit and an accurate model of the Earth (the WGS-84 ellipsoid). 
The results were compared to other results published, indicating a different final position of 
the aircraft from the locations of the published search area. 
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1. INTRODUCTION. On 8 March 2014, the flight Malaysian 370 from Kuala 
Lumpur to Beijing disappeared without any trace. The aircraft was a Boeing 777- 
200ER (777-2H6ER, ICAO type designator B772) with the registration 9M-MRO. 
After 32 minutes of flight preparation, the aircraft took off and flew normally for 49 
minutes. After the crucial point of IGARI at the border between Malaysian and 
Vietnamese airspace, the aircraft made abnormal manoeuvres for 1 hour and 17 
minutes (turns, climbs, descents) then settled for a long range flight of 5 hours and 
42 minutes. The total time flown abnormally was 7 hours, most probably until the 
fuel was exhausted. The long range en route flight of almost 6 hours was probably 
under automated guidance, very probably following a keyboard input flight plan 
changed earlier during the flight. Although all communication systems were shut 
down in a time window of 12 minutes (from 1709Z to 1721Z), the SATCOM satellite 
transponder continued to be functional and provided a number of very important clues 
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2 OCTAVIAN THOR PLETERAND OTHERS VOL. 69 
Table 1. MH370 Flight Phases. 
From To Phase Duration No 
2014 March 8 00:32:13 (1632Z) Flight preparation 32’ 0" 1 
00:00:13 (1600Z) 
E  ---------------------- Start of flight 32’ 
00:32:14 (1632Z) 00:40:38 (1641Z) ` Push-back, start and taxi 8' 26" 2 
00:40:38 (16412) — 01:01:16 (1701Z) Take-off and climb 20' 38" 3 
01:01:16 (1701Z) 01:21:04 (1721Z) En route to IGARI 19' 48" 4 
EE En End of normal flight | Start of abnormal 49’ 
manoeuvers under military surveillance 
01:21:04 (1721Z) 02:22:00 (1822Z) | Manoeuvers under military surveillance 1h0' 56” 5 
E En End of abnormal manoeuvers | Start of unknown 1h 1* 
flight 
02:22:00 (1822Z) 02:38:00 (1838Z) Unknown continuation of the flight, including Max 16’ 6 
? a turn onto a new flight track 
02:38:00 (1838Z) 08:19:29 (0019Z) Unknown en route flight Min 5h 7 
? 42' 
08:19:29 (0019Z) — 08:28:29 (0028Z) Glide descent with both engines out due to fuel  z 9" 8 
? ? starvation, until crash 
EE En End of flight 6h 07’ 


about flight MH370. This paper is the result of the original work of the authors, cal- 
culating all trajectories of the aircraft consistent with a single turn after loss of ground 
radar contact. Also, an optimisation cycle was used to calculate the most probable tra- 
jectories, which fit best on all the data known. The results are surprising, because they 
indicate a position where the flight could have ended, which is different from the past 
and the current search areas. Besides its scientific relevance, the authors hope that their 
work will prove useful in the future search operations. 


2. MH370 FLIGHT FACTS. Malaysian Airlines 370 took off with 239 Persons On 
Board (POB) (227 passengers out of the capacity of 282 and 12 members of crew). The 
initial flight plan was from Kuala Lumpur International (KLIA) to Beijing Capital 
International Airport (ZBAA). The ETA (Estimated Time of Arrival at planned des- 
tination (UTC)) was 2230Z and ETE (Estimated Time En route (estimated duration of 
flight in hours and minutes)) was 05:49. The fuel load was 49,100 kg, enough for over 
eight hours of flight. The MH370 flight phases are described in Table 1. Table 2 pro- 
vides a detailed sequence of known facts about the flight, including the seven relevant 
ping exchanges between the INMARSAT-3 Indian Ocean Region (IOR) satellite and 
the SATCOM transceiver, which allow the calculation of the flight trajectory. The key 
moments of the flight are defined in Table 3. 

The signal delay is called BTO (Burst Timing Offset) and it is caused by the time the 
radio waves take to travel from the satellite ground station to the satellite to the aircraft 
and back. Based on this, the instantaneous distance between the satellite and the air- 
craft (range) may be calculated for each ping exchange. The calculations presented in 
this paper rely mainly on these BTOs. 

The locus of the points within a certain range is a sphere with the centre in the sat- 
ellite position, with the range as the radius. This sphere may be intersected with the 
surface parallel to the ellipsoid or the geoid, at an altitude above the geoid, which 
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NO.1 RECONSTRUCTING THE MALAYSIAN 370 FLIGHT TRAJECTORY 3 
Table 2. Sequence of Known Facts about MH370. 

Time Event of the Flight MH370 Reference 

2014 March 8 INMARSAT I-3 IOR Log-On Avionics Power Up 


00:00:13 (16007) 
00:07:00 (16072) 


00:25:53 (16262) 
00:27:27 (16272) 


00:29:00 (1629Z) 
00:32:13 (16322) 
00:36:30 (16372) 
00:40:38 (16412) 


00:41:43 (16422) 
00:42:05 (16422) 
00:42:07 (16422) 
00:42:52 (16432) 


00:46:51 (16472) 
00:46:51 (16472) 
00:50:06 (1650Z) 
00:55:00 (16552) 
01:01:16 (1701Z) 
01:07:49 (17082) 
01:07:55 (170872) 
01:09:00 (1709Z) 
01:19:24 (17192) 
01:21:04 (17212) 
01:21:13 (17212) 
01:22:00 (17222) 
01:30:00 (17302) 


01:37:00 (17372) 


02:03:42 (18042) 
02:22:00 (18222) 
02:25:28 (18252) 
02:25:34 (18252) 
02:27:04 (18272) 


02:28:06 (1828Z) 
02:28:15 (1828Z) 


First Air DATA-2 ACARS! message BFO 86-90 Hz 
BTO 14780-14900 us 

“Delivery Malaysian Three Seven Zero good morning” 

“Ground Malaysian Three Seven Zero good morning 
Charlie One requesting push and start" 

ACARS message BFO 98-100 Hz BTO 14920-14940 us 

*Malaysian Three Seven Zero request taxi" 

"Tower Malaysian Three Seven Zero morning" 

TWR: “Three Seven Zero Three Two Right Cleared for 
take-off good night" MH370: “Three Two Right Cleared 
for take-off Malaysian Three Seven Zero thank you bye” 

Take-Off RWY 32R KLIA 

“Departure Malaysian Three Seven Zero” 

Cleared to climb to FL180 Direct IGARI 

Transferred from LUMPUR APP to LUMPUR 
RADAR 

“Lumpur Control Malaysian Three Seven Zero” 

Cleared to climb to FL250 

Cleared to climb to FL350 

ACARS message BFO 155-156 Hz BTO 15200-15220 us 

Pilot reports maintaining FL350 (unnecessary, 
unsolicited) 

HDG 025 GS 471 kts 

Last Air DATA-2 ACARS message BFO 131 Hz BTO 
15620 us 

Pilot reports maintaining FL350 (unnecessary, unsolicit- 
ed, repeated) 

FMS flight plan changed from the keyboard 

Transferred from LUMPUR RADAR? to HCM ATCC? 
on 120-9 MHz - Last transmission from cockpit: 
“Good Night Malaysian Three Seven Zero" 

Fly over IGARI N06?55'15" E103?34'43" (as tracked by 
LUMPUR RADAR) FOR) for 8 hours of flight 

Transponder switched off or failed 

Turn Left HDG 

HCM ACC requests other aircraft to relay communica- 
tions to MHA 370 — no answer 

Expected ACARS message not transmitted 

Climb to FL450 

Descend to FL120 (or lower) 

No response to Air DATA-2 ACARS - ACARS failed or 
has been switched off between 01:07:49 and 02:03:41 

200 NM NW of Pinang - Last Military Primary Radar 
contact 

Ping Log-On request from aircraft BFO 142 Hz BTO 
17120 us 

Ping handshake BFO 273 Hz BTO 51700 us 

Ping handshake BFO 175-176 Hz BTO 12560-12520 us 

Ping handshake BFO 144 Hz BTO 12500 us 

Ping handshake BFO 143 Hz BTO 12480 us 


First ACARS 


ATC Delivery contact 
LUMPUR GND contact 


LUMPUR TWR contact 
T/O clearance 


LUMPUR APP contact 


LUMPUR ACC contact 


Strange human action 


Last ACARS 
Strange human action 


(unconfirmed) 
Last COM 


Exit Malaysian airspace 


Last XPDR? 
PSRÓ data 
No COM 


No ACARS 

PSR data 

PSR data 

Proof of no ACARS 


Last PSR contact 
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4 OCTAVIAN THOR PLETER AND OTHERS VOL. 69 
Table 2. (Continued) 


Time Event of the Flight MH370 Reference 

02:39:53 (1840Z) Ground initiated Air Telephony call 

02:40:56 (1841Z) No answer — Ping handshake BFO 88-90 Hz Strange human action 
02:41:00 (1841Z) ACC Subang notifies about losing the flight Aircraft on a new track 


03:41:03 (1941Z) Ping handshake BFO 111 Hz BTO 11500 us 

04:41:05 (2041Z) Ping handshake BFO 141 Hz BTO 11740 us 

05:30:00 (21302) Search and Rescue (SAR) operation activated by Kuala 
Lumpur Rescue Coordination Centre (KL RCC) 

05:41:27 (2141Z) Ping handshake BFO 168 Hz BTO 12780 us 

06:41:22 (2241Z) Ping handshake BFO 204 Hz BTO 14540 us 


07:13:58 (2314Z) Ground initiated Air Telephony call Last Telephony Call’ 

07:15:02 (2315Z) No answer - Ping handshake BFO 216-219 Hz 

08:11:00 (0011Z) Ping handshake BFO 252 Hz BTO 18040 us Last complete Ping 

08:19:29 (0019Z) Last incomplete ping requested by user BFO 182 Hz Last incomplete Ping / 
BTO 23000 us both engines out 


08:30:00 (0030Z) Flight time limit based on FOB at FL310 or above at 
normal cruise airspeed TAS = 470 kts 
09:10:00 (0110Z) Expected Ping handshake did not arrive Proof of avionics out 


! ACARS: Aircraft Communications Addressing and Reporting System, automated bidirectional system 
for technical and operational messages between aircraft and a ground communication network or a satellite 
communication network 

? ACC Subang, Malaysia 

3 ACC Ho Chi Mingh, Vietnam 

^ FOB: Fuel on Board 

5 XPRD: Secondary Surveillance Transponder 

€ PSR: Primary Surveillance Radar (military) 

7 Not answered, but the subsequent ping proved aircraft to be in operational status. 


Table 3. Key moments of the MH370 flight. 


Malaysian Time UTC Identifier Description 

00:00:13 1600Z tL Avionics power up and satellite Log on 

00:41:43 1642Z tr Take-Off RWY 32R KLIA 

01:07:49 1708Z ta Last Air DATA-2 ACARS message 

01:19:24 1719Z tc Last COM “Good Night Malaysian Three Seven Zero” 

01:21:04 1721Z tx Last XPDR fly over IGARI N06°55'15” E103?34'43" 

02:22:00 1822Z tM Last Military PSR contact 200 NM NW of Pinang 

02:25:28 1825Z ty 1* complete satellite ping considered in the literature 

tu + TTT tm + TTT ts Start turn on a new track 

03:41:03 1941Z to 2"4 complete satellite ping, the first on the new track 

04:41:05 2041Z 13 34 complete satellite ping, the 24 on the new track 

05:41:27 2141Z t4 4" complete satellite ping, the third on the new track 

06:41:22 2241Z ts 5™ complete satellite ping, the fourth on the new track 

08:11:00 0011Z t6 eh complete satellite ping, the fifth on the new track 

08:19:29 0019Z t; 7 fg 7 incomplete user request ping — End of powered flight, 
start glide with both engines out 

tg + TGE to + TGE te End of flight 


corresponds with the presumed flight level of the aircraft. This is in fact the isobaric 
surface for level flight, corresponding to the local meteorological data. The intersec- 
tion yields a curve very close to a circle. At the given moment of the ping, the position 
of the aircraft is somewhere on that quasi-circle. 
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The shift in carrier frequency is called BFO (Burst Frequency Offset). It is caused by 
the Doppler-Fizeau effect due to the relative speed between the aircraft and the satel- 
lite. In the calculations below, the BFO values were not used, and this is a major dif- 
ference between this method and other methods used until now and known to the 
authors (Ashton et al., 2014; Hradecky, 2014). The reason for avoiding the use of 
BFO values is that the problem of relative speed is geometrically ill conditioned in 
this case, especially in the first half of the flight. 


3. CALCULATION OF QUASI-CIRCLES. This paper models the Earth as the 
WGS-84 ellipsoid. The ranges define quasi-circles at the intersection between the iso- 
baric surface that the aircraft levels on, and the sphere with the centre in the satellite 
antenna position and the BTO-based range as the radius. 

The INMARSAT-3 satellites are five satellites in 36,000 km high geostationary orbit 
(NO°), as follows: AOR-W at a longitude of W54°, AOR-E at W15-5?, I-3 at 25?E, IOR 
(the satellite used by MH370) at E64-5?, and POR at E178°. INMARSAT-3 IOR has 
an orbit inclined at 1:66?. This inclination gives a wobble movement with a period of T 
= 24 hours. The actual orbital altitude of IOR at ie was: 35,794 km, as compared to the 
target of 35,786 km Above Mean Sea Level (AMSL). Although geostationary, the 
actual position of the IOR satellite moved slightly during the MH370 flight. 

In the Earth-Centred-Earth-Fixed (ECEF) frame, the system of equations used to 
calculate the seven quasi-circles (see Figure 8) are the WGS-84 rotation ellipsoid equa- 
tions and the satellite-centred sphere equations: 


6378136-6 
Nog = (1) 


1 — e» x (sing)? 


x = (n, + H) cos p cosà 
y = (ny + H) cosgsinA (2) 
Z= Dell elt Hsing 


V = xsar) +(y ES ysar)? +(z = zsar) un Range satellite to aircraft — 0 (3) 
From (Ashton et al., 2014): 


1 : 
Range satellite to aircraft — 5€ E (BTO E bias) SE Rangesarellite to Perih GES (4) 
bias = —495679 us (5) 


The Equations (2) in x, y, z are substituted in Equation (3), providing a way to calculate 
longitude 4 as a function of latitude g. The seven quasi-circles correspond to the seven 
pings, illustrated in Figure 8, which occurred at the times t; to t; (see Table 3). 
However, the first complete satellite ping considered in the literature was not used in 
this paper, because the initialisation of the problem is uncertain, so less weight should 
be given to the starting point. Consequently, the work presented in this paper relies on 
the last six quasi-circles only. As a confirmation of this selection of quasi-circles, the con- 
tribution of the six pings to the total inconsistency of the trajectory is well balanced. The 
individual contribution of the first ping to the total inconsistency would have been dis- 
proportionately larger than the contributions of the other pings. 
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4. SOLVING TRAJECTORY EQUATIONS. The authors independently calcu- 
lated the MH370 trajectory based on the following data: 


1. Flight trajectory until the last known position P_18_22 as recorded by the mili- 
tary primary radars (N 6° 36’ 15” E 96° 33’ 14”) (Malaysian Gov., 2014; 
Hradecky, 2014). 

2. The set of recorded INMARSAT (INMARSAT, 2014) satellite pings (signal 
handshake delays and carrier frequency drift) (Malaysian Gov., 2014, 
Hradecky, 2014). 

3. The wind vector field, isobaric surface height and air temperature in the wider 
area at the moment of the flight, retrieved from the National Oceanic and 
Atmospheric Administration (NOAA) 2014). Interpolations in time, position, 
and flight levels were used for the best local current estimate of the wind velocity 
and the air temperature. 

4. The Boeing 777-200ER flight performance and systems (Boeing, 1992a). A 
Boeing 777-200ER flight simulator based on a dynamic model was developed 
by the authors (Section 5). 

5. The Boeing 777-200ER FMS (Flight Management System) specifications (Bulfer 
and Gifford, 1999). 

6. The WGS-84 ellipsoid shape approximation of the Earth. 

7. The INMARSAT-3 IOR satellite orbital data (Malaysian Gov., 2014). 


The unknowns of this problem and the space of solutions are considered in the fol- 
lowing subsections (see Figure 1). 

4.1. The Time to Turn Unknown (TTT). It was assumed that the flight continued 
straight for a while (on the track of 291° true) from the moment it vanished from the 
military primary radars (tm or P. 18 22). The duration of this straight segment TTT 
(Time To Turn) is one of the unknowns of the problem. TTT search domain was ini- 
tially considered 0-50 minutes, then reduced to 5-16-75 minutes, based on the 
following: 


e In numerical search trials, for TTT over 17 minutes, the breakdown of the quad- 
ratic error (£) on the six quasi-circles starts to show an abnormal error increase for 
the first two quasi-circles; 

e In both the uncorrected (Hradecky, 2014), and the corrected BFO diagram 
(Ashton et al., 2014), the turn is marked by a significant drop in the frequency 
due to the Doppler-Fizeau effect between the moments 1820Z and 1830Z, i.e. 
in an 11 minutes time interval; 

e The dynamic model simulation presented further provides an additional valid- 
ation for the upper limitation of this time to 16:75 minutes; if the flight continued 
on the initial track for 17 minutes or more, the optimised trajectory solutions (TK, 
FL and M) would largely disagree with the dynamic simulation (the fuel would 
have been exhausted before the last complete ping due to the required thrust set- 
tings resulting from FL and M); 

e No consistent solution emerged for TTT shorter than 5 minutes. 


After TTT, a coordinated turn was simulated, using a dynamic model Boeing 777- 
200ER simulator (see Section 5). The search step was one minute, refined down to five 
seconds for the promising domains. 
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Figure 1. MH370 trajectory problem unknowns: TTT; TK, FL, and M. 


This paper is based on the assumption that a single turn was made at the beginning 
of a long constant track segment. The BFO analysis (Figure 9 from Ashton et al., 2014) 
is consistent with a major turn between certain time moments: 1820Z and1830Z. After 
that turn, the BFO diagram indicates a linear variation of the Doppler-Fizeau fre- 
quency shift, i.e. no other turns. A notable fact is that one 360° turn at high M 
would typically last for 11-5 minutes, and the time window in the diagram is just 
too small for multiple turns concentrated at the beginning of the long constant track 
segment. On the other hand, multiple turns sparsely distributed over the entire 
voyage seem impossible given the same BFO diagram, because at each ping there is 
approximately the same BFO, the timing between pings is not always predictable 
(there are phone calls from the ground which spoil the one hour default time interval), 
and this is consistent with the hypothesis of an approximately constant ground speed 
vector. Another reason for the single turn assumption is the simplicity and the timing 

of the FMS programming. The programming of each additional turn into the FMS 
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should have been based on ad hoc navigation points entered individually as LAT- 
LONG. Attempting such a scenario on a flight simulator indicated a significant dur- 
ation of the operation, not consistent with the existing time window in this case. 

One additional rationale of the single turn hypothesis comes from the speculation on 
the goal of the person who did program the FMS as such. The goal of this action could 
have been to hide evidence in a place as remote as possible. The choice of track coin- 
cides with such a strategy, to end the flight as far as possible from any populated or 
frequented area. Programming multiple turns into the FMS would have been pointless 
and detrimental to the goal of flying as far as possible. This is speculative, but it only 
supplements all the other non-speculative grounds of the single turn assumption. 

4.2. The Track Unknown (TK). After TTT minutes, the aircraft took a turn to an 
unknown track TK (presumed constant), and from that moment it kept that track for 
the entire duration of the flight, until both engines became fuel starved. This is consist- 
ent with the variation of the frequency shift due to the Doppler-Fizeau effect on the 
carrier (the seven BFOs). Initially, all possible tracks were considered, (0?—359?), but 
a refined search was performed in the promising range of 170?—200? (the Southern 
range), and finally in the range 183°-193°, because all solutions obviously proved to 
be concentrated within these limits. The search step was 1°. The next best promising 
interval lays in the Northern range (around 336°) was not considered in the refined 
search for the following reasons: 


e The Northern tracks lead to inconsistencies at least three times larger than the 
Southern tracks; 

e In both the uncorrected (Hradecky, 2014), and the corrected BFO diagram 
(Ashton et al., 2014), the Southern variant matches very well the expected set 
of BFOs, in contrast with the Northern variant. 


This constant track phase of the flight (Phase 7 in Table 1) was probably flown under 
automated guidance, until the airplane became fuel starved, glided down and crashed. 
Other authors considered a constant heading for the flight. Although the Boeing 777 
autopilot has a function to maintain a constant heading (HDG HOLD on the Mode 
Control Panel-Boeing, 1992b), it is unlikely that this flight was not under FMS guid- 
ance (LNAV mode), and in this case HDG changes automatically to compensate for 
the crosswind. Thus, it is logical to assume a constant track flight (loxodrome or 
rhumb line) instead of a constant heading flight. A third hypothesis concerns a 
“DIRECT TO” a distant point flight, which would force navigation along the great 
circle (orthodrome). However, in the particular case of MH370, the difference 
between the orthodrome and the loxodrome routes is insignificant, and thus the con- 
stant track route covers both scenarios. A fourth hypothesis is a flight with the auto- 
pilot engaged in the “TRK HOLD” mode (constant track). This is also covered by 
the assumption made in this paper. 

4.3. The Flight Level Unknown (FL). The flight level (FL) of the cruise flight is 
presumed constant. Initially, all possible flight levels from 0 ft to 43,100 ft were consid- 
ered. The problem is not very sensitive to the geometric FL variable, except for the way 
the speed was formalised (see the Mach number below), where the temperature matters 
to a large extent. Also, the wind velocity changes from one level to the other, leading to 
significant differences between the flight levels. Very soon in the initial phase of the 
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Table 4. The subspace of solutions for refined search. 


Unknown Significance Units Min Max Step Layers 
TIT Time to turn Minutes 5 16-75 0-25 48 
TK Track i 183 193 1 11 
FL Flight level x100 feet 340 430 10 10 
M Mach number — 0-82 0:89 0-01 8 
Total 42,240 


optimisation process, it became obvious that flight levels below FL340 were impos- 
sible, for two reasons: 


e The inconsistency of these solutions is unacceptably high; 
e The dynamic model simulation indicated that achieving a time of flight until the 
last half ping (t7) would require a lower density of air. 


The fine phase of the optimisation was performed on a search range between FL340 
(34,000 ft) and FL430 (43,000 ft), with a step of 1,000 ft. 

The authors found no grounds to consider a variable flight level and Mach number 
along the cruise. Flying under FMS guidance would require a constant flight level 
selected for cruise, and the same value input in the altitude selector on the Mode 
Control Panel (Boeing, 1992b; Bulfer and Gifford, 1999). 

4.4. The Mach Number Unknown ( M). The Mach number (M) is presumed con- 
stant. The initial search domain was between M0:65 and MOO. but the envelope was 
correlated with FL following the flight performance of the aircraft type. The reason for 
choosing Mach instead of True Airspeed (TAS), Indicated Airspeed (TAS), or Ground 
Speed (GS) is related to the operational logic of the Flight Management System, Auto- 
Pilot and Auto-Throttle System on board the Boeing 777 aircraft. There is a slight dif- 
ference in the final result between maintaining a constant M and a constant TAS, and 
this comes from the outside temperature variations. The height of the isobaric surfaces, 
air temperatures and the wind velocity vectors along the route were taken from the 
NOAA archives (NOAA, 2014), as per the date and the hour of the flight. 
Interpolation in time, interpolation in the vertical profile of the wind vector field 
and interpolation in position data provide the best approximations of the real meteoro- 
logical parameters. For the later stages of the optimisation, the Mach search range was 
reduced to M0-:82-MO0:89. Whereas the maximum Mach number was delivering many 
consistent solutions, the lower speeds proved to be constantly inconsistent and were 
discarded. The search step for the Mach number was 0:01. 

4.5. Assumptions for the Last Segment of the Flight. The last incomplete user 
requested ping at t; is considered to mark the end of the powered flight, and the 
start of a powerless glide. The accurate simulation of the flight using a dynamic 
model (see Section 5) confirms this hypothesis, with a low error margin between the 
timing in the simulator and the timing of the ping signals (see Table 7). 

The last segment of the flight was a glide with both engines out due to fuel starvation 
(Pleter and Constantinescu, 2014). 

4.6. Trajectory Equations. Given the above considerations, the space of solutions 
was reduced from 34,056,000 initial possible solutions down to 42,240 possible trajec- 
tories, as summarised in Table 4 and illustrated in Figure 1 by underlined text. 
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Figure 2. Local wind triangle solved at each time step in the course of flight simulation. 


Considering the True North (u) and the True East (v) components of the local hori- 
zontal wind velocity vector retrieved from the NOAA records, the speed triangle 
system of trigonometric equations is the following (see Figure 2): 


TAS -cosHDG+v=GS -cos TK (6) 
TAS -sin HDG + u = GS - sin TK 


The solution of the system is the following: 


+TAS - cos TK + TAS? — (u- cos TK — v- sin TK)? 


HDG = ¥2tan | - 
Sn u-cos TK + (TAS —v)- sin TK (7) 


GS = v- cos TK + sin TK F \/ TAS? — (u-cos TK — v: sin TK)? (8) 


The system of six equations that result from the data (the time delays BTO of the six 
relevant pings) is over determined, and it was solved so as to minimise the quadratic 
sum of errors, which was considered as the inconsistency of the solution (E). The incon- 
sistency was calculated and used as an optimisation criterion to refine the search to the 


most plausible solution: 
e= Xd} — min (9) 


Each of the six errors d; is the distance between the aircraft at the moment of each of 
the six pings, and the closest point belonging to the respective quasi-circle. 

The flight was simulated using a Runge-Kutta IV method with a time step of 
10 seconds, integrating the local speed triangle (Figure 2) and the dynamic flight 
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Figure 3. The flight profile of the MH370 under military primary radar surveillance. 


simulator equations. At the moments of the pings, the position resulting from the simu- 
lation was compared to the position resulting from the radio transmission by calculat- 
ing the distance d; between the two. The quadratic sum of distances d; was minimised 
using a complete search optimisation method. All calculations were done by a 
MATLAB application written by the authors. 

The wind velocity and air temperature were retrieved from the NOAA world 
weather database (NOAA, 2014), with a granulation of 0-5? of latitude and longi- 
tude. Linear interpolation of latitude and longitude was used between the nodes. 
Vertically, the data are presented for the following isobaric surfaces: 10, 20, 30, 
50, 70, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 
800, 850, 900, 925, 950, 975, and 1000 mbar. These isobars are converted into 
flight levels using the ISA model, and linear interpolation was used between the 
isobaric surfaces. 


5. THE DYNAMIC MODEL FLIGHT SIMULATOR. For the purpose of this 
research, a Boeing 777-200ER flight simulator was developed. The simulator is 
based on a dynamic model i.e. considers the mass of the aircraft and calculates the 
fuel flow at each time step. The simulator served three separate goals: 


e To integrate the flight trajectory, the current mass of the aircraft, and the fuel 
consumption. By comparing the fuel consumption to the known fuel on board 
in the first stage of the flight, one can determine the approximate moment of 
fuel exhaust. This is compared to ig, which results from optimisation (see 
Table 7). 

e To calculate the fuel consumption of a dynamic climb to FL450, a descent to 
FL120 or below, and a climb back to an unknown cruise FL (Figure 3). 
These manoeuvres were recorded by the military primary surveillance radars, 
so the exact levels or details of the manoeuvres are not known. Instead of a 
simulator based on a time step integration, the authors used a method based 
on energy and aircraft climb / descent efficiency to calculate the fuel consumed 
by it. 

e To simulate the final glide descent portion, since the start of glide descent (ig) 
until crash (tg). 
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Table 5. Aerodynamic Coefficients for Boeing 777 (Pereira, 2011). 


Cia Coo Riet Kz 


54 0:013 —2-8° 0:047 


The equations of the simulator are the following: 


Cr = Cig: (4 — acro) (11) 
Cp = Cpo + K2- C (12) 
L=®. TAS! cO (13) 
D 
= — 15 
cos(a + T) MS) 
T-sin(a-- T) + L =m- g - Cosi (16) 


The coefficients used for Boeing 777 are as presented in Table 5. 
Simulation of the fuel consumption between the last ACARS message and last mili- 
tary primary radar contact was done with the following equations: 


FFavg : (tu — ta) + FCag = FC (17) 
AE, + AE, 
F Eu MAD iki 1 
CAE mm (18) 


The efficiency eff of the engines was calibrated based on the FOB at engine start and 
FOB at the moment of the last ACARS message (A). The taxi fuel consumption was 
calculated using the following formula: 


= Ataxi + Biazi ` Lat + Ctaxi : Na (19) 
T, amb 
where djaxi = —0:1223 braxi = 010335 Craxi = 0:1385 are the coefficients for a Boeing 777, 
na is the number of accelerations during taxi (it is considered to be at least two) and 
the ambient temperature in Kuala Lumpur in March during night was assumed 
to be Tamb = 297-15 K (Khadilkar, 2011). The taxi fuel consumption which results 
from the above equation is f=457:8 kg. The remaining FOB at take-off is 
49100-457-8 = 48612 kg. Knowing the FOB broadcast by the last ACARS 
message (43,800 kg) and considering an average fuel flow of FFavg = 6406 kg/h and 
the specific energy of Jet A-1, s= 43-02 MJlkg, the resultant efficiency of the engines 
is 32-32%. In this case, both the kinematic and the potential energy of the aircraft 
change. Using the calibrated engine efficiency, the resulting fuel consumption 
between the last ACARS message (A) and last military primary surveillance radar 
contact (M) is approximately 10,000 kg, depending on the FL to which the aircraft 
is considered to have climbed for the long cruise flight towards the south. For 
example, if the aircraft climbs to FL390, the calculated fuel consumption for the 
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Table 6. Equation (20) Coefficients. 


Dua P0 Pol 


31-47 0:004923 0:01895 


last ACARS - last primary radar contact segment is FC;,,;; = 10065-3 kg. In this case 
only the change in potential energy is considered (calculated with an average mass) and 
only at the climb regions. It is supposed that the fuel economy during the descending 
portions of the flight is balanced by the higher consumption in the dense air. 

Simulation of fuel flow during the long cruise flight is based on Boeing (1992a) and 
Filippone (2006). Using the calculated thrust and the mass of the aircraft, the follow- 
ing polynomial equation was defined (coefficients are given in Table 6): 


FF/ENG = poo + pio: mt poi: T (20) 


Simulation of the final glide descent with both engines out in the time interval TGD 
presumed 9 minutes is presented in Pleter and Constantinescu (2014). 

During the cruise flight, TK is kept constant by small turns, which compensate for 
crosswind. However, in the glide descent region, the simulation keeps a constant HDG 
instead of T'K, and these turns are no longer performed, as illustrated in the inset in 
Figure 8. This is consistent with the expected behaviour of the automated systems. 
By keeping a constant HDG, the local crosswind slightly drifted the aircraft during 
the gliding descent. 


6. ANALYSIS OF RESULTS. The multimodal optimisation using e as a criterion 
(inconsistency to be minimised) yields the most probable trajectories and the most 
probable values of the unknowns TTT, TK, and M. The unknown FL is considered 
a parameter, because the problem is not sensitive with respect to the geometric FL. 
Just two T'K ranges came out to be consistent: the Southern track around 190? 
(true) and the Northern track of around 336?. The Southern track fits better (a very 
low inconsistency, 38 solutions between 16 and 25 km), whereas the Northern track 
does not fit so well (zero solutions under 30 km). 

The multimodal optimisation yields the set of all solutions with a certain minimum 
probability instead of a single most probable solution. Given the nature of the problem 
and the scarcity of data coming from a unique source, this approach adds more value 
than the initial search done by the same authors for the most probable trajectory 
(Pleter and Constantinescu, 2014). Thus all trajectories with an overall inconsistency 
lower than 25 km were considered as potential solutions, and they are represented in 
Tables 7 and 8, and in Figures 4 to 6. The threshold of 25 km is remarkably low as com- 
pared to the total length of the trajectory (almost 6,000 km). There are 38 consistent 
solutions (with £ < 25 km), and they are analysed in Table 7 and in the subsequent 
Figures. In Figure 4 all the 38 consistent solutions were represented as cylinders. 
The diameter of the cylinder is the inverse function of the inconsistency, i.e. consistency 
(1/2). 

The most consistent solutions presented in Table 7 have a remarkably low inconsist- 
ency, of e < 17 km. Column “Fuel error” in Table 7 (second part) compares the top of 
descent (tg) calculated with the optimisation method based on the satellite pings and 
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Table 7. Consistent MH370 trajectory solutions and the corresponding positions of the end of flight (E) — solutions in yellow are confirmed by the available fuel criterion, 
with a fuel error < 1%. 


Solution & TIT TK M FL LAT, LONG, GS LAT; LONG; GS; LAT; LONG, GS; LAT; LONG; GS; 

(km) (mi) C) O Ou e C) (kt) ` C) IW (kt) (°) ©) (kt) (°) C) (kts) 
1 24:9 7:10 191 0:87 38,000 -2:58 93:51 5019 -10-79 91:88 508:1 —19-13 90:22 4935  —27-00 88-48 469-8 
2 23:3 7:25 191 0:86 36,000  —2:59 93-49 503:3 | —10:80 91:86 5065 -19-14 90:19 499-6  —27.02 88:46 470:2 
3 22:3 7:40 191 0-85 34,000 -2:61 93-47 5051 —10-81 91-83 5051 -19-16 90:17 5059 27-06 88-43 4711 
; n "s n "n 39,000 | -224 93-54 4983  —10:42 92-06 5052 -18:75 90:56 491:7  —26:61 88-99 467:5 
7 17:9 8:50 190 0:86 37,000 —2-27 93-52 499-2  —10:45 92-05 5037  -18 77 90:55 495.4 | —26:64 88-98 467-6 
8 16-4 9:05 190 0-85 35,000 —2-28 93-50 500:4  —10:46 92-02 5017  —18-78 90:52 500:9 26-67 88-95 468-4 
9 21:5 9:30 189 0-88 43,000 -1:92 93-63 489-6 9-95 92:31 4950 18-14 90:97 494-7 | —26:06 89-59 467:8 
13 162 10:25 189 0-86 38,000 —1-91 93-53 495.9  —10:06 92:21 50L3  -18:36 90:87 4938 26-23 89-47 465:4 
14 149 10:35 189 0-85 36,000 —1-904 93-52 496:3  —10:09 92:20 499-0  —18:38 90:85 497:3  —2626 89-45 466:4 
15 13:9 10:50 189 084 34,000 —1-96 93-49 497.2 —10-12 92-17 4969 — 18-40 90-82 501-1  -2630 89-42 467:8 
16 22:0 10:55 188 087 42,000 -1:64 93-66 487-4 —9:65 92:50 4924  —-1T82 91:31 491:9 —25-70 90-09 465:6 
17 247 11:00 189 0-86 37,0000 -1-88 93-45 499-2  —10:08 92:12 503:3  -1842 90-77 498:0 | —26:34 89-36 468-7 
18 19-7 11:10 188 0:87 42,000 —1-59 93-63 487-6 —9:61 92:47 4926 —17-77 91:29 492.0 | —25:66 90-06 466-0 
19 245 11:30 188 088 43,000 -155 93-57 490-2 —9:61 92-40 4948 —17-81 91:21 4953  -257] 89-08 469-6 
23 161 12:10 188 0-84 35,000 —1-63 93-5] 493-3 —9-74 92:34 494.2. —17:99 91-15 497.6 | —25:88 89-91 465:3 
24 23.5 12:30 188 086 38,000 -1:53 93-46 496-1 —9-70 92:28 501:3  -1802 91:09 4967 —25-92 89-84 466-7 
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25 18:0 12:40 187 0-87 43,000 -1:25 93-63 484-4 —9:22 92-62 488:1  -17:33 91:59 489:2 | —2524 90:53 467-9 
26 23.6 12:40 188 0-85 36,000 —1-55 93-43 496-2 —9-72 92-26 498:5  —18:03 91:07 499-] —25-95 89-82 467-6 
27 242 12:45 187 0:86 41,000 —1-28 93-63 485:5 —9:27 92-62 490:5  -1742 91:60 489-7 —25-28 90:53 464:1 
28 171 13435 187 087 42,000 -121 93-56 488-2 —9:24 92:54 492.5  -1742 91:50 4924 25-35 90-43 468:5 
30 21:5 13:35 187 085 38,000 -1:26 93-54 490-1 —9:34 92:52 495-6 | —17.58 91:49 493:0 —25:43 90:41 462:9 
31 21-6 13:40 187 0:84 36,000 —1-30 93-52 489-9 —9:38 92-51 492:5 -17-60 91:48 493:7 —25-47 90:40 463:3 
32 22-0 13:50 187 083 34000 -1:34 93-50 490-2 —9:43 92-49 489-7 —17:63 91:45 494-9 25-53 90:38 464:3 
33 23:9 14:00 187 0:86 39,000 —1-20 93-47 493-2 —9:33 92-45 499-4 -17-61 91-41 495-4  -2552 90:32 465:9 
34 20.7 14:45 186 087 443,000 —086 93-55 485-0 —8:85 92-68 4876 | —1697 91-80 489:3 —24-91 90-89 470:2 
35 194 14:55 186 086 41,000 —0-89 93-54 485-9 —8:90 92-68 490:4 —17-05 91-80 490066 —24-95 90-88 466:1 
37 23.6 15:45 185 0:86 43,000 —0-59 93-62 478-9 —8:49 92-90 480.7  —16:51 92:17 484:5 —24.39 91-43 467:2 
38 225 16:20 185 0:86 42,000 -0:56 93-54 482-4 —8-52 92-82 484:8 -16-61 92-09 487:8 | —2451 91:34 467-4 
Solution € TTT Tk M FL LAT; LONG; GS; LAT; LONG; GS; LAT; LONG; Fuel Error eh wlo wind 
(km) (min) © C mm © (s O O (s © © (%) ` Ge) (km 

1 24-9 7:10 191 0:87 38,000 -38:57 85-64 478-4 —39-66 85:34 475-1  -4042 85-14 —2:6% —681 195:8 
2 23:3 7:25 191 0:86 36,000  —38:54 85-63 479-0  —39-64 85:34 4740 —40-40 85-14 -46%  —1186 174-3 
3 22:3 7:40 191 0:85 34,000 -38:55 85-62 480.0 | —39-64 85:34 473:2 —40-41 85:15 -72% ——1859 154-6 
6 19-8 8:40 190 0:87 39,000  —38:17 86-41 4768 | —39:26 86.14 4736 | —40-02 85:96 -1:5% —392 194-7 
7 17-9 8:50 190 0:86 37,000 —38:14 86-41 4765 —39-24 86:15 4728 | —40:00 85:97 —3:2% —836 178-4 
8 16-4 9:05 190 085 35,000  -38-14 86-40 4770 39-23 86:14 4727] | —40:00 85:97 —5:57$  —1418 163-5 
9 21:5 9:30 189 0-88 43,000 —37-70 87-28 4768 —38-80 87.04 4743 —39-56 86:88 1:695 424 190-1 
13 162 10:25 189 0:86 38,000 —37-71 87.16 4753  -3881 86:93 471-1 —39-57 86:76 —2:1% —533 181-0 
14 149 10:35 189 0-85 36000 -37:71 87.16 4747  —38:80 86:93 4708 —39-57 86-77 —4:0%  — 1028 164-9 

Continued 
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Table 7. (Continued) 


Solution e TTT TK M FL LAT; LONG; GS; LAT; LONG; GS; LAT; LONG, Fuel Error E“ wlo wind 
km) (mm © C (d © 0O kis) © © kts) © © (4) ` Ge) (km 
15 139 10:50 189 0:84 34,000 —37-72 87:14 4746 | —38.82 86:92 470:9 | —39-58 86-76 —6:579 | —1688 144-7 
16 220 10:55 188 0:87 42,000 -37:25 88-05 4739 | —38:35 87:84 471:1 201 87-70 1:4% 355 197-6 
17 247 11:00 188 0:86 37,000 —37-88 87-05 477:2  —3899 86:81 4733 | —39-75 86:65 —3-2% —839 184-5 
18 19-7 11:10 188 0:87 42,000 -37:21 88-03 474] | —38331 87:82 471:3 | —39:07 87-68 1:4% 355 182-6 
19 245 11:30 188 0:88 43,000 —37-44 87-92 477:5 | —38:55 87-71 4750 | —3931 87-58 1:7% 427 204-7 
2 i 12:05 m n 37,000  —3727 87-91 4732  —38:36 87-70 469-0 | —39:13 87-56 —2:6% ES 170-4 
23 16:1 12:10 188 0:84 35,000  —3728 87-89 4723 38-37 87:69 467:9 | —3913 87-55 —4:9%  —1256 155-8 
24 235 12:30 188 086 38,000 —37-44 87:79 4761  —38:54 87:58 4721 | —3931 87-44 —2:1% —533 189-6 
25 180 | 12:40 187 0:87 43,000  —36:80 88-76 4734  —37.90 88:58 470:5 | —38:66 88-46 2:0% 517 184-7 
26 23:6 12:40 188 0:85 36,000 —37-44 87-78 4751  —38:54 87:58 471:0 | —39:30 87-43 -40% — 1034 173-8 
27 242 12:45 187 086 41,000 -3673 88-76 471:6 | —3782 88:58 468:4  —38.59 88-46 1:0% 268 1810 
28 171 13-15 187 087 42,000 -3693 88:65 4748 | —3803 88:47 4720 28 80 88-35 1:4% 358 194-7 
" m 13:35 " m" 38,000 | —36:83 88:64 471-8 | -37:92 88-46 467.9 | —38:69 88-33 -1:5% E. 174:3 
31 21-6 13:40 187 0:84 36,000  —36:84 88-64 470066 —37-93 88-46 4658  -38:69 88-33 —3-4% —869 158-3 
32 220 13:50 187 0:83 34,000  -36-87 88-63 469-7] —37-96 88-45 4640 | —3872 88:32 —5:99$  —1514 144-8 
33 239 14:00 187 0:86 39,000 -37-01 88:54 474.9 | —38:12 88:36 471:5 | —38:88 88-23 -10% —254 192-4 
34 207 14:45 186 (0087 43,000 —36-51 89-37 4746 | —3T.61 89-22 471:3 | —3837 89-12 2:0% 520 196-5 
E B BH e ees 41,000 —36-44 ps e —37:53 e E —38:30 Eum 1:0% SK Ba 
37 236 15:45 185 0:86 43,000 -35:87 90:18 471:1 — —36:97 90-05 4667  —3773 89-97 2:3% 604 185-9 
38 225 1620 185 0:86 42,000 -36:01 90-09 4721 -37-11 89-96 46719 | —37:88 89-88 1:8% 461 186-8 
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Table 8. Flight level distribution of the consistent MH370 trajectory solutions. 


FL 340 350 360 370 380 390 400 410 420 430 
Number of solution £ < 25 km 3 2 4 3 4 4 3 4 6 5 
Average € 19:41 16:22 20-87 19-86 21-54 21-32 18-73 20-53 21-57 21-65 
e 3 22 8 1636 2 2331 7 1794 1 24-93 6 19-83 5 2079 10 21-49 4 2433 9 21-46 
15 13-92 23 1608 14 1493 17 2466 13 1624 21 1699 12 1764 20 17:05 11 2388 19 2449 
32 2201 26 2364 22 1698 24 2345 33 2388 29 1775 27 2418 16 2203 25 18-01 
31 21-61 30 21-54 36 24-58 35 1941 18 1965 34 20-67 
28 1708 37 23-63 
38 2247 
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Figure 4. Position of the end of flight (£) for each of the 38 consistent solutions. 
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Figure 5. The locus of the end of flight positions (£*) for solution 77 considering a constant angular 
velocity turn during glide descent; the circles represent the possible points of impact for each 
respective angular velocity. 
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Figure 6. The end of flight point (F) for each of the 38 consistent solutions (discs); the weighted 
average of these points is represented by a light colour disc; the transparent areas resulting from 
the turning glide descent model (Figure 5) were positioned on the calculated end of flight points 
to give a probable area of impact; the ellipses on the 7'^ arc represent the most probable areas of 
impact. 


the top of descent (ig) calculated with the integration of the fuel flow until the 
moment of fuel exhaust. The relative error ranges between —7:2% and 2:39». For 
nine solutions (represented in yellow), the relative time error stays within +1% (254 
seconds of flight in absolute value), which is remarkable, given the inherent uncertainty 
caused by the manoeuvres before the last turn. Moreover, these small relative errors are 
symmetrically distributed in the range of the consistent solutions. Such a distribution 
appears to validate the results of the optimisation. The relative time error filter discards 
the lower FL solutions, such as solutions 3, 75, and 32. Also, this filter discards the 
lower M solutions, such as 8 and 26. Given the double validation, from the optimisa- 
tion and the fuel burn simulation, solutions at the upper end of the FL and M scales 
stand out, e.g. solution 71. 

The places where the flight might have ended in the 38 most probable scenarios are 
represented in Figures 4 and 5. The places seem clustered in seven isolated groups, but 
the gaps between these groups might be equally interesting for searching. The solution- 
clustering phenomenon is probably a consequence of selecting integer values for TK. In 
reality, TK could have very well been a rational number, in the case of the third hypoth- 
esis of the FMS programming presented above, a “DIRECT TO” a distant point given 
by LAT and LONG. Figure 8 illustrates the whole simulated trajectory, which corre- 
sponds to one of the solutions (17). The last portion of this trajectory 1s made 
visible at a better scale in the inset. 

These 38 end of flight points were calculated assuming that the glide descent was 
flown with a constant heading (no turn). Since the probability of a turn due to the 


Downloaded from https://www.cambridge.org/core. IP address: 191.96.87.184, on 09 Sep 2018 at 18:56:58, subject to the Cambridge Core terms of use, available at 
https://www.cambridge.org/core/terms. https://doi.org/10.1017/S0373463315000570 


20 OCTAVIAN THOR PLETERAND OTHERS VOL. 69 


$22°30 


T+opic-of-Capricorm. 


Google earth 


Figure 7. The end of flight point (£) for each of the 38 consistent solutions, compared to the current 
search areas and to other results published (CNN, 2014; Australian Government, 2014; Ashton 
et al., 2014; Pleter and Constantinescu, 2014). 


initial thrust asymmetry is high (ATSB, 2014a; 2014b - MH370 - Definition of 
Underwater Search Areas), and the rate of turn (the angular velocity) is not known, 
a locus of all potential end of flight points (E*) was developed and illustrated in 
Figure 5. All trajectories in Figure 5 correspond to the end of flight point of 
Solution 17, and are all turning trajectories with incremental angular velocities. The 
circles at the end of each same-length orthodrome trajectory represent potential 
points of impact with the water. 

In reality, the turn could have been with variable angular velocity. The fixed angular 
velocity model in Figure 5 however includes in the envelope of the E* points all pos- 
sible impact points. In Figure 6, for each cluster of E solutions, the corresponding E* 
envelope was represented as a transparent area. Out of this area, the ellipses on the 
seventh arch illustrate the spiral dive, which could be a higher probability impact area. 

In Table 8 the same 38 solutions were placed on their FLs, making 1t obvious that the 
fuel burn criterion favours the higher FLs solutions (in yellow). Flying in a denser at- 
mosphere, although geometrically consistent with the BTOs, would have fuel starved 
the engines much earlier than the moment tg, which is known. 

For the same solution 17, the simulated ground speed variation over the entire flight 
was represented in Figure 9. The ground speed is calculated as the vector addition 
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Figure 8. Quasi-circles of intersection between the satellite-centred spheres and the isobaric surface 
of the aircraft taken parallel to the WGS-84 rotation ellipsoid and the trajectory of MH370 as 
simulated for solution 77. Inset illustrates final segment with a slight drift both ways to the effect 
of the crosswind on the glide descent trajectory. 


between the horizontal component of the true airspeed and the local wind velocity. For 
a possible validation of such a trajectory with a method based on BFO, Table 9 indi- 
cates the ground speed in the seven moments of the seven pings, decomposed in the 
ECEF framework. 

The inclusion of the wind vector field in simulating the flight increased the relevance, 
as shown in Table 7. The inconsistency of the best-fit trajectories is between 6 and 15 
times worse in cases where the wind vector field is not used. This validates the simula- 
tion and the use of the proper wind vector field, which is one of the few certain facts 
regarding this flight. 
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Figure 9. Simulated MH370 ground speed variation (solution 17). 


Table 9. Solution 11 ground speed decomposed in ECEF reference for the seven quasi-circles. 


ECEF 1 2 3 4 3 6 7 


x'(m/s ` 24610194 41:08357 42:09721 41-21286 39-08501 33-95766 33-19933 
y’ (m/s) 12:56416 —5:61979  —42.02263  —TT-TTI39  —104-78226  —150:01248  —152-87285 
Z'(m/s) 113:88568  —248.75327  —248.00685  —238.06273 .—212-.53179  —190-78845  —186:75045 
t (s) 210 4743 8346 11967 15564 20940 21450 
gs (kts) 528:3 493-4 498-9 496:3 469-3 478:3 475-4 


7. CONCLUSIONS. The method presented uses the ranges determined by the last 
six pings along the unknown portion of the cruise flight, avoiding calculations based on 
the carrier frequency drift. Thus, an issue of geometric ill conditioning of the problem is 
avoided (since the relative aircraft-satellite speed vector is almost perpendicular to the 
aircraft ground speed vector). In a multimodal optimisation loop with complete 
search, the flight is simulated based on a dynamic model, determining the unknowns 
of the flight, which best fit the known facts. Thus, the method gives simultaneously 
the time from the last known position until the left turn, the track, and the Mach 
number of the cruise. The flight level can be used as a parameter of the optimisation. 

All the 38 solutions found have a flight ending point E falling closely to each other, as 
illustrated in Figures 4, 6, and 7. The search area proposed by this paper may be approxi- 
mated as a rectangle aligned with the 7th arch, with the following corners: S39?46'35" 
E84?35'51"; S40?43'18" E85?24'09"; $37?38'09" E91?09'40"; S36?47'14" E90?12'56". 

Figures 6 and 7 illustrate the most probable end of flight points calculated in this 
paper, in contrast to other published results (Ashton et al., 2014; Hradecky, 2014; 
CNN, 2014; Australian Government, 2014; Pleter and Constantinescu, 2014). 
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